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摘 要 ;估计 寨 卡 和 登革热 (DEN) 疫 情 从 亚洲 输入 中 国 及 引起 本 地 暴发 的 风险 。 基 于 2015—2017 年 国外 疫情 和 流动 人 口 
数据 ,构建 输入 模型 估计 输入 病例 数 ,并 计算 基于 不 同 温度 和 群体 免疫 水 平 下 分 支 过 程 的 本 地 疫情 传播 概率 及 基本 再 生 
数 。 中 国 的 寨 卡 输入 病例 主要 来 自 新 加 坡 、 泰 国 和 越南 ,预测 的 病例 数 分 别 为 7.0C(95% CI:6. 5 一 7.5)、2.0(95%% CI:1.8 一 
2.2) 和 1.0(95% CI:0.9 一 1. 1) ;登革热 输入 病例 主要 来 自 泰国 、 马 来 西亚 新加坡、 越南 、 非 律 宾 、 印 度 尼 西亚 、 印 度 和 韩 
了 亲 , 预 测 的 病例 数 分 别 为 700. 0(95% CI:679. 8~720. 2) .654. 1(95% CI:641. 8 一 666. 2)、376. 3(95% CI: 368. 2 一 384. 1), 
27. 1(95% CI:268. 55~285. 33) 241. 2(95% CI:233. 6 一 248. 8) .67. 0(95% CI:59. 6~74. 5),9. 1(95% CI:6.7~11. 3) 和 
-830(95% CI:1.9 一 4.1)。 温 度 在 28. 9 左右 是 最 适宜 寨 卡 和 登革热 传播 的 条 件 , 此 时 发 生 本 地 传播 的 风险 概率 分 别 为 
4% 和 99.9%。 将 人 类 群体 免疫 水 平 从 0 增加 到 0. 2 和 0. 6, 寨 卡 和 登革热 的 基本 再 生 数 分 别 为 8.1、6.7、3.2 和 3.2、 
用 7.1.3。 输 入 病例 更 多 来 自 于 南亚 ,中 国 中 南 和 东南 是 暴发 本 地 传播 的 高 风险 地 区 ,特别 是 在 6-8 月 。 新 加 坡 疫情 更 容易 
导致 寨 卡 在 中 国 传播 ,而 泰国 .越南 、 马 来 西亚 和 新 加 坡 疫 情 是 中 国 发 生 登 革 热 本 地 传播 的 最 大 导 火 索 。 
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Abstract; The aim of this paper is to estimate the risk of Zika and dengue fever (DEN) infection that imports from Asia and 
causes local outbreaks. Based on the oversea epidemic data and the floating population between 2015 and 2017, an input 
model was built to estimate the number of imported cases, and the local epidemic transmission probability and basic regener- 
ation number were calculated via the branching process under different temperatures and herd immunity levels, The impor- 
ted cases of Zika in China mainly came from Singapore, Thailand and Vietnam, with the predicted number of cases as 7. 0 
(95% CI:6. 5-7.5), 2.0 (95% CI: 1. 8-2. 2) and 1.0 (95% CI:0. 9-1. 1). The imported DEN cases were mainly from Thai- 
land, Malaysia, Singapore, Vietnam, Philippines, Indonesia, India and Korea, the predicted number of cases as 700. 0 
(95%CI;679. 8-720. 2), 654.1 (95% CI: 641. 8-666. 2), 376.3 (95% CI; 368. 2-384. 1), 277.1 (95% CI; 268. 6-285. 3), 
241. 20 (95% CI:233. 6-248. 8), 67.03 (95% CI:59. 6-74. 5), 9.1 (95% CI: 6. 7-11. 3) and 3.0 (95% CI:1. 9-4. 1). The 
optimum temperature for Zika and DEN transmission is around 28. 9 °C , in which the risk probability of local transmission is 


24.4% and 99.9%, respectively. When the human herd immunity level is 0, 0. 2 and 0.6, the basic reproduction numbers 
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of Zika and DEN are 8.1, 6.7, 3.2 and 3.2, 2.7, 1.3, respectively. The imported cases mainly come from South Asia. 


South-central and south-eastern China are top-risk areas for local transmission, especially in June-August. The infection in 


Singapore is more likely to cause Zika outbreak in China, while the infection in Thailand, Vietnam, Malaysia and Singapore 


are the biggest cause of local transmission of DEN in China. 


Key words; Zika; dengue fever; risk assessment; imported cases; local transmission 


近 几 十 年 来 ,人 类 旅行 的 大 幅 增长 和 范围 扩大 导 
致 传染 病 在 全 球 蔓 延 ,典型 的 例子 包括 蚊 媒 传染 病 和 
最 近 的 新 冠 肺炎 5 。 寨 卡 和 登革热 主要 通过 伊 蚊 传 
播 ,集中 分 布 在 热带 和 亚热带 地 区 ,亚洲 是 主要 疫 区 
之 一 沾 。2015 一 2016 年 , 寨 卡 病毒 暴发 了 前 所 未 有 
的 爆炸 性 疫情 ,席卷 南美 ,并 草 延 至 亚洲 5 ,世界 卫 
生 组 织 记录 显示 ,共有 87 个 国家 和 地 区 报告 本 地 案 
卡 病毒 感染 证 据 ,疑似 和 确诊 病例 超过 140 万 中 。 近 
几 年 登革热 疫情 暴发 的 频率 和 规模 在 持续 上 升 ,全 球 
约 2 弛 8 个 国家 面临 风险 ,每 年 大 概 3. 90 亿 人 被 感染 ， 
由 9 600 万 例 出 现 临床 症状 ,造成 极为 严重 的 疾病 
PA), 2019 年 世界 卫生 组 织 将 登革热 列 为 全 球 
健康 威胁 之 一 ""。 对 这 两 类 传染 病 在 中 国 本 十 
揪 的 风险 评估 具有 重要 的 科学 意义 。 
= 据 国家 旅游 局 、 中 国旅 游 研究 院 等 有 关机 构 数 据 
#42015 年 中 国 出 境 旅游 人 数 达 到 1. 17 亿 人 次 ， 
旅游 人 数 达 到 1. 34 亿 人 次 ,其 中 亚洲 地 区 占 
68M0% 。Lai 等 "发 现 东南 亚 到 中 国 的 航空 旅客 数 
量 从 2005—2015 年 翻 两 番 , 随 着 航空 运输 量 增加 , 东 
南亚 登革热 传人 中 国 的 可 能 性 显著 增加 。 岳 玉 娟 
等 中 探讨 了 2014—2018 年 中 国 登革热 重点 疫 区 病 
例 罕 间 特 征 及 相关 因素 ,发 现 登革热 病例 集中 在 广东 
和 甘南 省 境内 ,本 地 病例 与 境内 输入 病例 ,境外 输出 
病 鲍 高 度 相关 。 刘 小 波 等 发 现 2019 年 中 国 绝 大 
多 数 项 目 省 媒介 伊 蚊 密度 超过 登革热 传播 或 暴发 阔 
值 ,登革热 本 地 传播 及 暴发 风险 高 。 杨 志俊 等 5 回 
顾 了 中 国 输入 病例 的 发 现 过 程 . 旅 行 史 以 及 确诊 途径 
等 流行 病 学 要 素 ,提出 需要 加 强人 境 人 口 的 得 查 , 并 
注意 控 蚊 的 建议 。Zhang 等 中 也 提出 要 及 时 监测 从 
疫情 发 生 国 入境 人 员 感 染 情况 的 建议 。 

目前 ,从 亚洲 地 区 输入 寨 卡 和 登革热 的 病例 数 与 
中 国 本 地 寨 卡 和 登革热 传播 程度 及 在 时 空 上 的 潜在 
风险 尚 不 清楚 。 鉴 于 此 ,整合 亚洲 疫情 、 人 口 旅游 情 
况 和 气候 环境 等 数据 ,借助 MATLAB 进行 数值 分 
析 ,评估 疫情 输入 中 国 和 引起 本 地 传播 的 风险 。 采 用 
Dorigatti 等 "中 开发 的 数学 模型 估计 2015—2017 年 
亚洲 地 区 流行 期 间 预计 引入 中 国 的 寨 卡 和 登革热 感 
染 人 数 。 温 度 在 蚊 媒 病毒 感染 中 起 关键 作用 ,其 影响 
包括 调节 病毒 和 蚊虫 的 传播 因子 (转化 率 、 蚊 虫 叮咬 


率 和 死亡 率 )。 量 化 寨 卡 病毒 和 登革热 输入 中 国 的 风 
险 ,并 使 用 Johansson 等 25 描述 的 分 支 过 程 方法 佑 
计 输 入 寨 卡 和 登革热 后 本 地 传播 概率 。 


1 材料 与 方法 


1.1 数据 来 源 


基于 中 国 知 网 、 百 度 学 术 、 必 应 学 术 和 SCE 
HUB 文献 检索 库 输入 登革热 (dengue fever), 38 F 
(zika) 和 风险 评估 (risk assessment) 等 关键 词 搜 寻 相 
关 文 献 ,按照 相关 度 进 行 筛选 。 从 世界 卫生 组 织 中 获 
取 2015 一 2017 年 寨 卡 和 登革热 的 病例 数 与 流行 时 间 
窗口 (https://ntdhq, shinyapps. io/dengue5/)。 中 
国 出 境 的 旅游 人 数 在 国家 统计 年 鉴 中 获得 Chttp:// 
www. stats. gov. cn/tjsj/ndsj/) , 新加坡、 菲律宾 、 韩 
E ERBE 印度尼西亚、 泰国 、 马 来 西亚 .越南 和 日 本 等 
地 区 的 入 境 旅 游人 数 从 中 国 国家 统计 局 -文化 和 旅游 
部 获得 (https://data. stats. gov. cn/gjwz. htm 和 ht- 
tps://www. stb. gov. sg/content/stb/en. html), 
2015—2017 年 全 国 698 个 气象 站 点 每 日 气象 资料 在 
中 国 科 学 数据 共享 服务 平台 获取 (http://data. cma. 
cn) ,将 各 省 份 所 含 站 点 的 月 平均 值 作为 各 省 月 气温 
数据 。 从 百度 获取 各 个 国家 2015 年 的 人 口 数 据 (ht- 


tp://www. hotelaah. com/ Asia. html) 。 
12 寨 卡 和 登革热 输入 病例 数 的 风险 分 析 


根据 Dorigatti 等 " 引 开 发 的 数学 模型 ,将 流行 病 
学 数据 与 旅行 数据 相 结 合 ,以 量化 疾病 传 入 风险 。 该 
模型 通过 入 境 和 出 境 模拟 一 种 疾病 传人 未 受 影响 国 
家 的 2 种 方式 。 设 O( 亚 洲 地 区 除 中 国 ) 代 表 疾 病 的 
起 源 地 ,S (中国 ) 代 表 未 受 疾病 感染 的 地 区 。 入 境 病 
例 数 是 指 由 S 地 前 往 O 地 旅行 的 游客 ,在 O 逗留 期 
间 被 感染 ,然后 在 恢复 之 前 返回 S 的 游客 人 数 。 故 


S BO 的 入境 病 例 数 
Ts+90 =T s+oAP y » (1) 
其 中 ， 
NoLo 
“二 OP 天， (2) 
Py =min( EL, ih (3) 
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出 境 病例 数 是 指 在 O 地 受 感染 ,在 潜伏 期 或 染 (a) an 
病 期 到 达 S 地 旅游 的 人 口 数量 。 故 O 到 S 的 出 境 病 ie 
例 数 其 中 ,0 过 ;过 1。 本 地 传输 的 复合 函数 为 
Bow = Nob str Y setel a 58s) 
p g6 =(1+ 54 g6) 
TW 其 一 个 特性 是 g, (0) 代表 零 代 灭绝 的 概率 。 因 此 ， 
Ps=365POP,’ (5) 输入 一 个 病例 和 多 个 病例 的 本 地 传播 概率 可 分 别 表 
=, [TetTi IRN 
Pi min wW 1). (6) Po et pC Os (14) 
即 输 入 中 国 的 病例 数 Pe =1— ge Os (15) 
N =Is:0 + Eoas (7) 若 人 类 种 群 由 于 先前 的 暴露 或 交叉 保护 而 并 非 
1. 3 案 卡 和 登革热 本 地 传播 的 风险 评估 模型 完全 匈 感 , 则 输入 病例 后 天 绝 的 可 能 性 会 增加 。 为 


此 ,使 用 有 效 繁殖 数 Rw 进行 调整 , 即 用 Rw 替换 
基本 再 生 数 Ro 表示 在 完全 易 感 人 群 中 由 一 个 ” Rw 。 参 数 RY 是 通过 调整 Rw 与 免疫 ( 原 免 疫 ) 人 群 
原 发 病例 产生 的 次 级 病例 的 平均 数量 ,是 传播 性 的 关 ”的 比例 而 得 出 De ， 


键 奖 量 标准 > OFA Johansson 等 "中 描述 的 分 支 过 R? =RY(1— P rm). (16) 
BRIER, SAE HED. 学 参数 依赖 温度 ,是 温度 T 的 函数 ea BGR 
= Ro=RoRy ， (8) ~ Herp, UR BE T 取 某 确定 值 时 , 式 (17)、 
© R? =GoBiw¥T 1» (9) (8)、(19) 和 (20) 的 值 随 之 确定 , 即 可 确定 模型 (9) 和 
N RY =a, ao 《10) 中 的 参数 值 。 
© u(T) =0.869 2—0. 159T +0. 0111 6T? 一 
se. Ro 表示 每 个 感染 者 产生 感染 载体 的 平均 数 0. 000 340 8T3 +0. 000 003 809T', (17) 
vuh fa. Fu > yu ` STA 
ENR 表示 每 个 传染 性 载体 产生 的 传染 性 人 类 的 平 a(T) 0.004 3T +0.0943, GD 
bi er 20 Bro (T) =max{0. 001 044T(T — 
an Ey Mail ap he 0.072 9T —0. 903 7 
Wiest MSH ko WORE ASS AY HERE A PRD Bu (T) | 12.4 C<LT<L26.1 C; (20) 
I g,(s) 和 gi(s) 来 描述 , 即 1, 26.1 C=T<32.5 C., 
K an hu =$ 
© = ap 


表 1 亚洲 多 国 与 寨 卡 和 登革热 疫情 相关 的 数据 


寨 卡 病 登革热 
地 区 名 称 人 口 数 量 
寨 卡 病 发 生年 病例 数 / 例 登革热 起 止 年 病例 数 / 例 

新 加 坡 5 607 283 2016 411 2015—2017 27 203 
印度 尼 西 亚 (印尼 ) 258 383 264 2015 1 2015 一 2017 373 227 
菲律宾 105 172 921 2017 57 2015—2017 543 350 
泰国 68 714 519 2015 686 2015 一 2017 262 073 
印度 1 338 658 835 2017 3 2015 一 2017 417 480 
日 本 126 994 512 2016 8 

韩国 51 217 803 2017 19 2015 一 2016 574 
马来西亚 (大 马 ) 31 660 700 2016 2 2015—2017 294 153 
越南 93 640 435 2016 219 2015 一 2016 207 610 


注 : 寨 卡 起 止 时 间 具 体 如 下 :2016 年 7 月 1 日 到 12 月 31 日 (新 加 坡 );2017 年 1 月 1 日 到 2 月 29 日 (菲律宾 );2017 年 1 月 1 日 到 2 月 29 日 
(印度 );2016 年 9 月 1 日 到 12 月 31 日 (马来西亚 ) ;一 表示 无 报告 病例 数 (日 本 ) 。 
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表 2 ”模型 参数 描述 


参数 定义 取 值 来 源 
POPo 亚洲 地 区 的 人 口 总 数 v 表 1 
Ts 每 年 亚洲 地 区 到 中 国 的 旅行 者 人 数 v 国家 统计 局 .中 华人 民 共 和 国文 化 和 旅游 部 
T so 疫情 期 间 中 国 到 亚洲 地 区 的 旅行 者 人 数 v 国家 统计 年 鉴 
No 在 亚洲 地 区 的 病例 数 / 例 v #1 
Ww 流行 时 间 窗 口 /天 v 表 1 
A 在 亚洲 有 逗留 的 人 均 感染 风险 式 (2) 文献 [15] 
Py 在 潜伏 期 或 染病 期 返回 亚洲 其 他 国家 的 概率 /% 式 (3) 文献 [15] 
Ps 在 疫情 窗口 期 间 亚 洲 其 他 国家 到 中 国 的 人 均 概率 / % 式 (5) 文献 [15] 
Pi 感染 者 在 流行 时 间 窗 口 处 于 潜伏 期 或 染病 期 的 概率 /% 式 (6) 文献 L15] 
1/4 平均 病 媒 寿命 /天 式 (17) 文献 [19-20] 
a 每 只 蚊子 每 天 叮咬 的 次 数 式 (18) 文献 [19] 
genic 人 到 病 媒 的 有 效 传播 率 式 (19) 文献 [19-21] 
Pu 病 媒 对 人 的 有 效 传播 率 式 (20) 文献 [19-21] 
od 旅行 者 在 亚洲 地 区 的 平均 逗留 时 间 / 天 7 文献 [22] 
Cr 感染 寨 卡 的 人 类 内 在 潜伏 期 /天 3. 50 文献 [23] 
Cr， 感染 寨 卡 的 人 类 染病 期 /天 15 文献 [24] 
S.. 感染 登革热 的 人 类 内 在 潜伏 期 /天 5 文献 [25] 
= 感染 寨 卡 的 人 类 染病 期 /天 6 文献 [25] 
C, 人 均 具 性 蚊子 数量 /只 0. 35 文献 [26] 
= 病 媒 生物 存活 外 部 潜伏 期 的 平均 比例 /% 0. 89 文献 [27] 
Ni 色散 参数 0. 10 文献 [28] 


Le 模型 验证 


(6 基于 模型 估计 从 9 个 受 寨 卡 和 登革热 影响 的 亚 
洲 国家 (新 加 坡 、 非 律 宾 , 韩 国 ,印度 .印度 尼 西亚 、 泰 
国 S 名 来 西亚 ,越南 和 日 本 ) 疫 情 波及 到 中 国 的 风险 。 
模 骨 验 证 在 MATLAB R2018a 运行 ,结果 均 通 过 检 
验 ,对 应 的 了 小 于 0.05。 假 设 游客 在 海外 去 留 的 平 
均 天 数 (7 天 ) 代 表 来 自 亚洲 地 区 的 游客 在 中 国 的 平 
BEA? 。 假 设 负 二 项 概率 生成 函数 ,色散 参 
数 固定 在 二 0. 1 处 ,以 估计 本 地 传输 的 概率 55。 候 
设 有 症状 的 寨 卡 和 登革热 病例 与 健康 人 有 相同 的 旅 
行 概率 ,这 9 个 国家 与 中 国 的 旅行 者 人 数 通 过 国家 统 
计 年 鉴 与 各 个 国家 的 旅游 局 获得 。 汇 总 9 个 国家 疫 
情 时 间 窗口 (从 2015 年 以 后 第 一 批 确诊 病例 出 现 症 
状 到 2017 年 最 后 一 批 确诊 病例 之 间 的 天 数 ) 内 的 旅 
RIM. 

为 验证 模型 的 有 效 性 ,进行 案例 研究 。 输 入 病例 
数 预期 数量 的 不 确定 性 是 通过 从 潜伏 期 和 感染 期 分 
布 中 抽样 1 000 个 来 计算 95% BMA Kia). HAE 
群体 的 免疫 力 对 疫情 起 重要 作用 ,探讨 当 人 类 免疫 能 


力 不 断 增强 时 , 寨 卡 和 登革热 在 中 国 的 传播 概率 。 利 
用 温度 对 蚊 媒 传染 病 的 影响 ,探讨 不 同 温度 和 免疫 条 
件 对 基本 再 生 数 Ro 的 影响 。 


2 结果 


图 1 为 亚洲 地 区 寨 卡 和 登革热 的 流行 导致 中 国 
的 潜在 病例 。 预 测 寨 卡 病毒 在 新 加 坡 输入 形成 中 国 
的 潜在 病例 (输入 和 输出 ) 数 量 最 多 ,为 7.00 例 (95% 
CI:6. 50~7. 50) ,其 次 是 泰国 2. 00 例 (95% CI:1.80 一 
2. 20) 和 越南 1. 00 例 (95% CI:0. 80 一 1. 20) , 相 比 之 
下 ,菲律宾 .印度 .二 国 、 马 来 西亚 ,印度尼西亚 和 日 本 
疫情 几乎 不 会 影响 到 中 国 ( 见 图 1(a))。 预 测 登 革 热 
在 泰国 输入 中 国 的 潜在 病例 数 最 多 700. 00 例 (95% 
CI:679.79 一 720.21), 其 次 是 马来西亚 654. 07 例 
(95% CI:641. 83 一 666. 24) 、 新 加 坡 376. 30 例 (95% 
CI: 368. 17 一 384.13)、 越 南 277.14 Pi) (95% CI: 
268. 55~285. 33) .菲律宾 241. 20 例 (95% CI:233. 61 
一 248. 79) 和 印度 尼 西 亚 67. 03 例 (95% CI: 59.57 一 
74. 45) ,而 印度 9. 05 例 (95% CI:6.71 一 11. 32) 和 韩国 
3.00 例 (95% CI:1. 90~4. 10)( 见 图 1(b))。 
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(b) 登革热 从 亚洲 输入 中 国 的 病例 数 


图 1 从 亚洲 地 区 输入 到 中 国 潜在 病例 数 的 估计 值 ,黑色 条 表示 95% 的 置信 区 间 


图 2 为 其 他 国家 寨 卡 和 登革热 疫情 导致 中 国 在 
不 同 群体 免疫 水 平 下 发 生 本 地 感染 的 概率 。 发 现 本 
地 暴发 的 概率 随 着 群体 免疫 水 平 的 增高 而 缓慢 降低 ， 
其 申 案 卡 的 降低 幅度 更 明显 ,但 均 需 要 群体 免疫 大 于 
890% 才 有 显著 下 降 。 当 群体 免疫 水 平 为 0 时 ,新 加 
坡 = 泰 国 和 越南 引起 中 国 发 生 寨 卡 本 地 传播 的 概率 分 
I 19.54% (95% CI:17.44% ~21. 64%) 4. 88% 
(95% CI: 4.60% ~ 5.16%) 和 1.29% (95% CI: 


149% ~1. 38%) ,而 菲律宾 、 印 度 、 韩 国 、 马 来 西亚 、 
ae 
So [ 一 一 新 加 a 
置信 区 间 = 
一 一 泰国 
a 置信 区 间 
=> h S 
425 f = = 
Se : Š 
w x 
K 
E 
rt i 
名 人 类 群体 免疫 力 
(a) ZF 


印度 尼 西 亚 和 日 本 的 寨 卡 本 地 病例 较 少 ,引起 中 国 本 
地 病例 的 概率 几乎 为 零 ( 图 2(a))。 与 寨 卡 相 比 , 4 
群体 免疫 水 平 为 0 时 ,菲律宾 印度尼西亚 .印度 和 昔 
国 引起 中 国 本 地 暴发 的 概率 分 别 为 93. 98% (95% 
CI: 89. 44% ~ 98. 53%) 59. 71% (95% CI:55. 31% ~ 
64. 10%). 9.86% (95% CI: 8.25% ~ 11. 48%) 和 
4.11% (95% CI:2. 87% ~5. 35%) GUA 2(b))。 而 
登革热 在 泰国 、 越 南 、 马 来 西亚 和 新 加 坡 的 疫情 导致 
中 国 本 地 暴发 的 概率 达到 了 99. 00%( 图 2(b))。 


80 | 度 尼 西 3 
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(b) 登革热 


2 不 同 群体 免疫 水 平 下 导致 本 地 传播 发 生 的 概率 ,阴影 区 域 表 示 95% 的 置信 区 间 


图 3 为 寨 卡 和 登革热 在 不 同 免疫 条 件 和 温度 下 
基本 再 生 数 R, 的 变化 。 发 现 随 着 温度 增加 , Ro 均 
先 增 后 减 , 呈 倒 U 型 。 温 度 在 28. 9 “CIN, BER AE 
EAW Ro 达到 最 大 值 ,分 别 为 8. 11 和 3. 24( 无 群体 
免疫 时 )。R。 随 着 群体 免疫 力 增加 而 逐渐 减 小 , 当 群 
体 免 疫 力 从 0 至 1 以 0.03 步 长 增加 时 ,车 群体 免疫 
力 从 0 增加 到 0. 20 和 0. 60( 见 图 3 箭头 处 ), 则 寨 卡 
和 登革热 的 基本 再 生 数 最 大 值 分 别 降低 为 6.65、 
3.24 和 2.66、1.29。 

图 4 为 寨 卡 和 登革热 在 不 同 温度 下 亚洲 各 国 疫 
情 对 中 国 本 地 传播 的 影响 大 小 。 结 果 表 明 ,当中 国平 
均 温 度 为 28.9 'C 时 ,新 加 坡 、 泰 国 、 越 南 和 韩国 的 疫 


情 导 致 中 国 发 生 寨 卡 本 地 传播 的 概率 分 别 为 
21. 14% (95% CI:19. 33%~22.94%)、5.34% (95% 
CI:4.53% ~ 6.14%), 1.47% (95% CI: 1.36% ~ 
1.58%) AN 0.39% (95% CI: 0. 36% ~0. 41%) (Al 4 
a). JEEE EN EE Pa, EN A eS BR 
生 登 革 热 本 地 传播 概率 分 别 为 95. 30% (95% CI: 
91. 03% ~ 99.57%), 60. 23% (95% CI: 54.39% ~ 
66. 06%). 10.56% (95% CI: 8.62% ~ 12.50%) 和 
3. 88% (95% CI:3. 23% ~4. 53%) (图 4(b))。 此 时 
以 上 国家 输入 中 国 的 总 病例 数 引起 本 地 传播 的 风险 
概率 分 别 为 24. 43% A 99. 92%。 当 中 国平 均 温度 下 
降 至 23.9 CC 时 ,以 上 国家 导致 中 国 发 生 寨 卡 本 地 传 
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图 3 不 同 温度 和 群体 免疫 力 下 寨 卡 和 登革热 的 基本 再 生 数 尺 , 的 变化 


播 概率 分 别 下 降 15.76%, 22.51%, 31.73% 和 37.77%。 当 中 国平 均 温 度 下 降 到 15.9 时 ,以 上 概 
15. 28%; 以 上 国家 导致 中 国 登 革 热 病 发 生 本 地 传播 ROP I) F E 85. 87%、85.01%、86. 42%, 87.47% All 
概率 分 别 下 降 14.87%, 16.23%, 28.13% 和 75. 41 88. 08% 79. 81% Fil 85. 18%, 


CO 
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© E 7=23.9°C : 
© MB 7=15.9°C es 
pa = 
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CN 新 加 波 泰国 越南 韩国 日 本 菲律宾 大 马 印度 “印尼 泰国 ”大 马 新 加 波 越南 菲律宾 印尼 印度 “韩国 
> (a) 案 卡 病毒 在 不 同 温度 下 的 本 地 传播 概率 (b) 登革热 在 不 同 温度 下 的 本 地 传播 概率 
x< 1 寨 卡 和 登革热 在 高 .中 和 低 风险 下 输入 到 中 国 的 本 地 传播 概率 ,本 地 传播 概率 显示 95% 置 信 区 间 


Sa 5 为 寨 卡 和 登革热 在 不 同 温度 条 件 下 在 中 国 高。 而 到 了 9 一 10 月 除了 广东 ,广西 和 海南 ,其 他 地 

市 的 传播 风险 ( 即 基 本 再 生 数 ) 。 结 果 发 现 ,每 年 区 的 基本 再 生 数 均 迅速 降低 。 在 11 月 到 次 年 3 月， 
6 一 8 月 是 疫情 的 高 发 期 ,其 中 广东 、 广 西 .海南 、 浙 ” 全 国 各 地 区 的 基本 再 生 数 均 小 于 1。 特 别 地 ,青海 和 
江 ,\ 福 建 , 江 西 和 湖南 等 南部 地 区 温度 较 高 , 风险 较 ” 西藏 地 区 全 年 的 基本 再 生 数 始终 为 0。 
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(a) 寨 卡 病毒 在 我 国 各 市 的 基本 再 生 数 O 登革热 在 我 国 各 市 的 基本 再 生 数 


5 我国 各 省 市 发 生 寨 卡 和 登革热 传播 的 逐 月 (4 一 11 月 ) 风 险 


3 讨论 传播 的 概率 具有 明显 的 差异 性 。 罕 卡 病毒 输入 中 国 
的 风险 主要 来 自 新 加 坡 ,其 次 是 泰国 和 越南 ,导致 本 
研究 发 现 不 同 国家 的 输入 风险 和 输入 导致 本 地 ”地 传播 的 概率 均 小 于 22%。 登 革 热 的 输入 风险 主要 
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来 自 泰国 越南 、 马 来 西亚 和 新 加 坡 ,其 次 是 菲律宾 。 
而 且 这 些 国家 的 输入 病例 很 大 程度 上 会 引起 登革热 
在 中 国 流行 , 概 率 均 超过 95%。 需 注意 的 是 中 国 至 
今 为 止 并 未 出 现 本 地 的 寨 卡 疫情 ,主要 原因 在 于 中 国 
部 署 了 严密 的 联防 联 控 策 略 以 应 对 2015—2017 年 亚 
洲 地 区 的 疫情 ”中 ,有 效 降低 了 寨 卡 病毒 在 中 国 地 
区 的 本 地 传播 风险 。 登 革 热 是 比 寨 卡 更 流行 的 蚊 媒 
传染 病 ,输入 中 国 病例 和 本 地 传播 风险 均 非 常 高 。 主 
要 原因 在 于 ,这 些 国家 的 登革热 病例 数 较 多 , 且 具 有 
大 规模 的 出 和 人 境 人 数 。 因 此 ,亚洲 地 区 暴发 的 寨 卡 和 
登革热 疫情 对 中 国 的 影响 不 容 忽视 。 

分 析 表 明 , 寨 卡 和 登革热 在 中 国 暴发 期 的 高 风 
险 时 间 是 6 一 8 月 ,集中 分 布 在 中 南 和 东南 部 ,在 11 
月 到 次 年 3 月 全 国 风 险 均 很 低 , 其 余 时 间 广 东 、 广 
西 和 海南 也 存在 一 定 风险 。 风 险 的 时 空 异 质 性 主 
要 时 温度 起 关键 的 调控 作用 ,其 对 基本 再 生 数 的 影 
mB ES 性 倒 UCU 型 ,最 适宜 温度 是 28.90 C. 5 
b PE (A HE fi 6 够 降低 输入 风险 ,但 只 有 当 群 体 免 
FE 60% 以 上 时 ,登革热 的 基本 再 生 数 才 有 可 能 小 
|。 实 际 上 ,登革热 在 中 国 是 因为 输入 病例 引起 
地 暴发 ,监控 数据 显示 ,每 年 的 2 一 7 月 均 有 不 


完 渍 例 从 南亚 和 东南 亚 进入 中 国 。 作 为 常住 人 口 
和 流动 人 口 最 多 的 广东 省 ,是 输入 风险 最 高 的 地 


DED 。 而 近 两 年 因为 对 新 冠 的 管控 作用 ,限制 了 人 
口 葛 大 范围 流动 ,相应 地 ,登革热 的 输入 和 本 地 感 
淆 后 人 数 大 幅度 降低 。 这 也 表明 人 口 流动 是 导致 
ere ae ee 


一 基于 输入 病例 及 出 人 Be AR. iH at a SC t E, 
MATLAB 数值 分 析 等 方法 ,评估 了 寨 卡 和 登革热 这 
2 类 收 媒 传染 病 通过 航空 旅行 从 亚洲 流行 国家 输入 
到 中 国 的 风险 以 及 引起 本 地 传播 的 概率 ,并 在 中 国 城 
市 的 尺度 上 估计 了 蚊 媒 的 时 空 分 布 。 本 研究 融合 了 
温度 ,人口 流动 和 国外 疫情 等 更 多 动态 因素 。 研 究 结 
果 可 为 评价 中 国 寨 卡 和 登 娃 热 的 流行 趋势 和 风险 评 
佑 提供 科学 依据 。 

提出 的 数学 模型 可 以 与 高 分 辩 率 向 量 丰 度 和 和 气 
候 适宜 性 估计 相 结合 ,研究 方法 可 应 用 于 分 析 其 他 新 
出 现 或 再 次 出 现 的 蚊 媒 传染 病 。 本 研究 主要 存在 以 
下 局 限 性 : 人 
模型 未 吉 括 所 有 潜在 因素 ,例如 降雨 量 、 社 会 经 济 和 
随机 因素 等 。 


王 而 营 等 , 寨 卡 和 登革热 输入 中 国 及 引起 本 地 传播 的 风险 信 计 。 IT 
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